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We present an update of results from the search for microlensing towards the Large Magellanic Cloud (lmc) 
by EROS (Experience de Recherche d'Objets Sombres). We have now monitored 25 million stars over three years. 
Because of the small number of observed microlensing candidates (four), our results are best presented as upper 
limits on the amount of dark compact objects in the halo of our Galaxy. We discuss critically the candidates and 
the possible location of the lenses, halo or LMC. We compare our results to those of the MACHO group. Finally, 
we combine these new results with those from our search towards the Small Magellanic Cloud as well as earlier 
ones from the EROSl phase of our survey. The combined data is sensitive to compact objects in the broad mass 
range 10"^ — 10 M©. The derived upper limit on the abundance of stellar mass machos rules out such objects as 
the dominant component of the Galactic halo if their mass is smaller than 2M0. 



1. Research context 

The search for gravitational microlensing in our 
Galaxy has been going on for a decade, follow- 
ing the proposal to use this effect as a probe of 
the dark matter content of the Galactic halo [Q . 
The first microlensing candidates were reported 
in 1993, towards the lmc and the Galactic 
Centre by the EROS, MACHO and OGLE collab- 
orations. 

Because they observed no microlensing candi- 
date with a duration shorter than 10 days, the 
EROSl and MACHO groups were able to exclude 
the possibility that more than 10% of the Galactic 
dark matter resides in planet-sized objects 

However a few events were detected with longer 
timescales. In their two-year analysis [Q, the 
MACHO group estimated, from 6-8 candidate 
events towards the LMC, an optical depth of or- 
der half that required to account for the dynami- 
cal mass of a "standard" spherical dark halcQ; the 
typical Einstein radius crossing time of the events, 
tE, implied an average mass of about 0.5 Mq 
for halo lenses |l^. Based on two candidates, 
EROSl set an upper limit on the halo mass frac- 
tion in objects of similar masses |ll],|7||, that is 
below that required to explain the rotation curve 
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of our Galaxy^. 

The second phase of the EROS programme was 
started in 1996, with a ten-fold increase over 
EROSl in the number of monitored stars in the 
Magellanic Clouds. The analysis of the first two 
years of data towards the Small Magellanic Cloud 
(SMC) allowed the observation of one microlensing 
event ||l^ also detected by This single event, 
out of 5.3 million monitored stars, allowed eros2 
to further constrain the halo composition, exclud- 
ing in particular that more than 50 % of the stan- 
dard dark halo is made up of 0.01 — O.SMq objects 
jl^ ]. In contrast, an optical detection of a halo 
white dwarf population was reported jist , com- 
patible with a galactic halo full of white dwarfs. 

Very recently, the MACHO group presented an 
analysis of 5.7 year light curves of 10.7 million 
stars in the lmc with an improved determi- 
nation of their detection efficiency and a better 
rejection of background supernova explosions be- 
hind the LMC. They now favour a galactic halo 
MACHO component of 20% in the form of 0.4 Mq 
objects. Within a few days, the detection of a 
halo white dwarf population at the level of a 10% 
component was also reported by |jl^. Simulta- 



^ Assuming the original two EROSl candidates are mi- 
crolensing events, they would correspond to an optical 
depth six times lower than that expected from a halo fully 
comprised of MACHOs. 
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neously, the eros2 group presented its results 
from a two-year survey of 17.5 million stars in 
the LMC [Q. One erosI microlensing candidate, 
EROSl-LMC-2, was seen to vary again, 8 years af- 
ter its first brightening, and was thus eliminated 
from the list of microlensing candidates. Two new 
candidates were identified (eros2-lmc-3 and 4). 
Because this is much lower than expected if MA- 
CHOS are a substantial component of the galac- 
tic halo, and because these two new candidates 
do not show excellent agreement with simple mi- 
crolensing light curves, EROS chose to combine 
these results with those from previous EROS anal- 
yses, and to quote an upper limit on the fraction 
of the galactic halo in the form of machos. 

In this article, we describe an update on the 
EROS2 LMC data, an analysis of the three-year 
light curves from 25.5 million stars. While the 
sensitivity is improved, the main conclusions are 
unchanged compared to iQ. One of the two- 
year candidates was seen to vary in the third sea- 
son and was thus rejected. Three new candidates 
have been detected. We combine these eros2 
LMC results with those of previous independent 
EROS analyses, and derive the strongest limit ob- 
tained thus far on the amount of stellar mass ob- 
jects in the Galactic halo. 



2. Experimental setup and LMC observa- 
tions 

The telescope, camera, telescope operation and 
data reduction are as described in 12|. Since 
August 1996, we have been monitoring 66 one 
square-degree fields in the LMC, simultaneously in 
two wide passbands. Of these, data prior to May 
1999 from 39 square-degrees spread over 64 fields 
have been analysed. In this period, two thirds of 
the fields were imaged about 210 times in aver- 
age; the remaining third were imaged only about 
110 times. The exposure times range from 3 min 
in the LMC center to 12 min on the periphery; 
the average sampling is once every 4 days (resp. 
8 days). 



3. LMC data analysis 

The analysis of the LMC data set was done us- 
ing a program independent from that used in the 
SMC study, with largely different selection crite- 
ria. The aim is to cross-validate both programs 
(as was already done in the analysis of erosI 
Schmidt photographic plates |ll|) and avoid los- 
ing rare microlensing eventsF. The analysis is 
very similar to that reported in [|l8|. We only 
give here a list of the various steps, as well as a 
short description of the differences with respect 
to our two-year analysis. A detailed description 
of the analysis will be provided in |^,|l[ . 

We first select the 6% "most variable" light 
curves, a sample much larger than the number of 
detectable variable stars. This subset of our data 
is "enriched" in genuine variable stars^, but also 
and mainly in photometrically biased light curves, 
i.e. those of stars especially sensitive to the ob- 
serving conditions, such as stars very close to neb- 
ulosities or to bright stars. Working from this 
"enriched" subset, we apply a first set of cuts to 
select, in each colour separately, the light curves 
that exhibit significant variations. We first iden- 
tify the baseline fiux in the light curve - basically 
the most probable fiux. We then search for runs 
along the light curve, i.e. groups of consecutive 
measurements that are all on the same side of the 
baseline fiux. We select light curves that either 
have an abnormally low number of runs over the 
whole light curve, or show one long run (at least 
5 valid measurements) that is very unlikely to be 
a statistical fluctuation. We then discard light 
curves with a low signal-to-noise ratio by requir- 
ing that the group of 5 most luminous consecutive 
measurements be significantly further from the 
baseline than the average spread of the measure- 
ments. We also check that the measurements in- 
side the most significant run show a smooth time 
variation. 

The second set of cuts compares the measure- 
ments with the best fit point-lens point-source 
constant speed microlensing light curve (hereafter 



We have checked that the present program finds the 
same SMC candidate as reported in flld]. 
^ We monitor our selection efficiency with Monte-Carlo 
simulated variable star and microlensing light curves. 
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"simple microlensing" ) . They allow us to re- 
ject variable stars whose light curves differ too 
much from simple microlensing, and are suffi- 
ciently loose not to reject light curves affected by 
blending, parallax or the finite size of the source, 
and most cases of multiple lenses or sources. We 
also require that the fitted time of maximum mag- 
nification lie within the observing period or very 
close to it, and that the fitted timescale is shorter 
than 300 days. The latter cut is equivalent to 
requiring that the baseline flux of the star is ob- 
served for at least a few months; this is necessary 
in any analysis using this baseline flux. At this 
stage of the analysis, all cuts have been applied 
independently in the two passbands. 

After this second set of cuts, stars selected 
separately in the two passbands represent about 
0.01% of the initial sample; almost all of them 
are found in two thinly populated zones of the 
colour-magnitude diagram. The third set of cuts 
deals with this physical background. The first 
zone contains stars brighter and much redder 
than those of the red clump; variable stars in 
this zone are rejected if they vary by less than 
a factor two or have a very poor fit to simple 
microlensing. The second zone is the top of the 
main sequence, where the selected stars, known 
as blue bumpers display variations that are 
almost always smaller than 60% of the base flux 
or at least 20% lower in the visible passband 
than in the red one. These cannot correspond to 
simple microlensing, which is achromatic; neither 
can they correspond to microlensing plus blend- 
ing with another unmagnified star, as it would 
imply blending by even bluer stars, which is very 
unlikely. We thus reject all candidates from the 
second zone exhibiting one of these two features 
(see Fig. |lD. 

Compared to the analysis in jl^ , two new cuts 
are introduced to reject other types of variable 
stars that were not present in the two-year anal- 
ysis. The first one is aimed at stars which have a 
roughly constant luminosity for some time, then 
vary typically over one or two months to reach 
a new constant level. We cannot yet conclude 
whether these are physical variable stars or some 
kind of instrumental problem. The second cut 
is aimed at novae and supernovae. It rejects light 
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Figure 1. A colour-magnitude diagram in one 
EROS field, showing the location of candidates 
identified as "blue bumpers" , either from their 
small amplitude, or from their chromatic varia- 
tion (larger in the red than in the visible pass- 
band). 



curves which have a rise time significantly smaller 
than the decline time; it is not applied to events 
with a timescale longer than 60 days, in order not 
to reject microlensing phenomena with parallax 
effects, that also show an asymmetry. 

The final cuts are simply tighter cuts on the fit 
quality, applied to both colours (whereas similar 
previous cuts were applied independently in each 
passband), and a requirement that the observed 
magnification be at least 1.40 . 

The tuning of each cut and the calculation 
of the microlensing detection efficiency are done 
with simulated simple microlensing light curves, 
as described in For the efficiency calcula- 

tion, microlensing parameters are drawn as fol- 
lows : time of maximum magnification tg uni- 
formly within the observing period ±150 days. 
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impact parameter normalised to the Einstein ra- 
dius uq G [0,2] uniformly, and timescale tE G 
[1,400] days uniformly in \n{tE)- All cuts on 
the data were also applied to the simulated light 
curves. 
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Figure 3. Same as Fig. 2 for candidates EROS- 
LMC-6 and 7. 



Figure 2. Light curves of candidates EROS-LMC- 
3 and 5 (visible passband) . The plain curves show 
the best point-lens point-source fits; time is in 
days since Jan 1.0, 1990 (JD 2,447,892.5). 



Only four candidates remain after all cuts. Of 
the two candidates presented in eros2-lmc- 
3 is still a member of this hst, while eros2-lmc-4 
was seen to vary at least twice in the third season 
and was thus rejected. There are three new can- 
didates, numbered 5 to 7. Their light curves are 
shown in Figs. |^ and 0; microlensing fit parame- 
ters are given in Table ]l|. Although the candidates 
pass all cuts, agreement with simple microlensing 
is not excellent. In particular, eros2-lmc-5 is 
dubious : it has a bad fit to simple microlensing 
and is located in an atypical region of the colour- 
magnitude diagram. The geometric mean of the 



candidates timescales is about 32 days, including 
that of the erosI candidate lmc-1. 

The microlensing detection efficiency of this 
analysis, normalised to events with an impact pa- 
rameter uq < 1 and to an observing period Tobs of 
three years, is summarised in Table |^. The main 
source of systematic error is the uncertainty in 
the influence of blending. Blending lowers the 
observed magnifications and timescales. While 
this decreases the efficiency for a given star, the 
effective number of monitored stars is increased 
so that there is partial compensation. This ef- 
fect was studied with synthetic images using mea- 
sured magnitude distributions ||2^ . Our final ef- 
ficiency is within 10% of the "naive" sampling 
efficiency. Compared to the efficiency in |l^ , the 
present one is improved for the longest and short- 
est durations, but slightly lower for average dura- 
tions around 50 days. This is largely explained by 
the fact that we have included in the present anal- 
ysis stars in external LMC fields that were sampled 
less frequently. 
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Table 1 

Results of microlensing fits to the four eros2 
LMC candidates; ts is the Einstein radius crossing 
time in days, uq the impact parameter, and c^{^' 
the R{V) blending coefficients, constrained to be 
smaller than unity. 
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Table 2 

EROS2 detection efficiency in % for the LMC 3- 
year analysis, as a function of the Einstein radius 
crossing time i^; in days, normalised to events 
generated with uq < 1, and to Tobs = 3 yrs. 



4. Limits on Galactic halo MACHOs 

EROS has observed six microlensing candidates 
towards the Magellanic Clouds, one from erosI 
and four from EROS2 towards the LMC, and one 
towards the SMC. As discussed in [|l2|, and fur- 
ther in p3|] , we consider that the long duration 
of the SMC candidate together with the absence 
of any detectable parallax, in our data as well as 
in that of the MACHO group |jl^, indicates that 
it is most likely due to a lens in the SMC^. For 
that reason, the limit derived below uses the five 
LMC candidates. (The limit corresponding to all 
six candidates would be about a factor 1.13 times 
the limit shown, for masses larger than 0.01 solar 
mass.) 

The limits on the contribution of dark com- 
pact objects to the Galactic halo are obtained 
by comparing the number and durations of mi- 
crolensing candidates with those expected from 
Galactic halo models. We use here the so-called 
"standard" halo model described in as model 
1, but have checked that we obtain similar results 
for other reasonable halo models. The model pre- 



® Alternatively, it can be argued that, if due to: 
halo lens, this event corresponds to a heavy lens |l3Ul4lt 



dictions are computed for each EROS data set in 
turn, taking into account the corresponding de- 
tection efficiencies (|11 |Jl^ and Table || above), 
and the four predictions are then summed. In this 
model, all dark objects have the same mass M; we 
have computed the model predictions for many 
trial masses M in turn, in the range [10~^ M0, 
102 Mq]. 

The method used to compute the limit is as 
in |ll). We consider two ranges of timescale Ie, 
within or outside the interval I = [7.5; 190] days. 
This interval was chosen as follows. We first de- 
termine the average mass corresponding to the 
mean duration of the five LMC candidates, at 
about 0.2 Mq. We then compute the expected 
distribution of microlensing timescales for this av- 
erage mass and check that the observed spread in 
timescales for the candidates is compatible with 
the width of this distribution. This means that 
our candidates are compatible with the hypothe- 
sis that their spread in mass contributes very little 
to the width of the timescale distribution. The in- 
terval I is then chosen as a symmetrical interval in 
In (is) that contains 99% of the timescale distri- 
bution for halo MACHOs of 0.2 M©. Of course, aU 
five LMC candidates have timescales well within 
the interval I. 

We can then compute, for each mass M and any 
halo fraction /, the combined Poisson probability 
for obtaining, in the four different EROS data sets 
taken as a whole, zero candidate outside I and 
five or less within I. For any value of M , the limit 
/max is the value of / for which this probability is 
5%. Whereas the actual limit depends somewhat 
on the precise choice of I, the difference (smaller 
than 5%) is noticeable only for masses around 
0.01 and 10 Mq. Furthermore, we consider our 
choice for I to be a conservative one. 

Figure^ shows the 95% C.L. exclusion limit de- 
rived from this analysis on the halo mass fraction, 
/, for any given dark object mass, M. The solid 
line corresponds to the five LMC candidates; it is 
the main result of this article. This limit rules 
out a standard spherical halo model fully com- 
prised of objects with any mass function inside 
the range [10~^ — 10] Mq. In the region of stellar 
mass objects, where this result improves most on 
previous ones, the new LMC data contribute about 
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Figure 4. 95% C.L. exclusion diagram on the 
halo mass fraction in the form of compact ob- 
jects of mass M, for the standard halo model 
(4x lO^^M© inside 50 kpc), from all LMC and SMC 
EROS data 1990-99. The soHd Hue is the limit in- 
ferred from the five LMC microlensing candidates. 
The new MACHO 95% C.L. accepted region is the 
closed contour, with the preferred value indicated 
by the cross [ p^ . 



73% to our total sensitivity (the SMC and EROSl 
LMC data contribute 10% and 17% respectively). 
The total sensitivity for t^; = 50 days, that is pro- 
portional to the sum of Tobs e(i_E — 50d) over 
the four EROS data sets, is about 3.2 times larger 
than that of ||l^] and two thirds that of p^ . 

5. Discussion 

After nine years of monitoring the Magellanic 
Clouds, EROS has a meager crop of five microlens- 
ing candidates towards the LMC and one towards 
the SMC, whereas about 30 events are expected 
towards the lmc for a spherical halo fully com- 
prised of 0.4 M© objects. Moreover, some of 



the candidates cannot be considered excellent. 
These candidates were obtained from four dif- 
ferent data sets analysed by independent, cross- 
validated programs. So, the small number of ob- 
served events is unlikely to be due to bad (and 
overestimated) detection efficiencies. 

This allows us to put strong constraints on the 
fraction of the halo made of objects in the range 
[10"'' Mq, 10 Mq], excluding in particular at the 
95% C.L. that more than 40% of the standard 
halo be made of objects with up to IM© . The pre- 
ferred value quoted in |l^], / = 0.2 and 0.4Mo, is 
consistent with our limit as can be seen in Fig. ^. 
(The upper part - about 25% - of the domain al- 
lowed by [ p^ is excluded by the limit we report 
here.) 

There are several differences which should be 
kept in mind while comparing the two experi- 
ments. First, EROS uses less crowded fields than 
MACHO with the result that blending is relatively 
unimportant for EROS. (Were EROS results to 
be corrected for blending, the detection efficiency 
would increase slightly and the reported limit 
would be stronger.) Second, EROS covers a larger 
solid angle (64 deg^ in the lmc and 10 deg^ in 
the SMC) than MACHO, which monitors primar- 
ily 15 deg^ in the central part of the lmc. The 
EROS rate should thus be less contaminated by 
self-lensing, i.e. microlensing of lmc stars by dim- 
mer LMC objects, which should be more common 
in the central regions. The importance of self- 
lensing was first stressed in [ p4|p^ . Third, the 
MACHO data have a more frequent time sampling. 

The results from EROS and MACHO are appar- 
ently consistent, but the way they are interpreted 
is different. MACHO reports a signal and consid- 
ers the contamination of its sample as low or null. 
EROS2 quotes an upper limit and does not claim 
its sample to be background-free. The position of 
the lenses along the line of sight, halo or Magel- 
lanic Clouds, is also an issue. MACHO has com- 
pared the spatial distribution of its candidates 
across the face of the LMC and observes a better 
agreement with the halo hypothesis than with a 
specific model of the lmc. On the other hand, 
because the EROS stars are spread over a wider 
field, the fact that the EROS sample corresponds 
to a lower central value of the event rate (about 



7 



twice lower than that of macho) is compatible 
with an interpretation where a noteable fraction 
of the events are due to self-lensing. The small 
number of EROS candidates precludes at present 
any definitive conclusion on that topic. 

It seems likely that the single most important 
input to the question of the position of the lenses 
will come from the comparison of the microlens 
candidates samples towards the SMC with those 
towards the lmc. Because the two lines of sight 
are rather close (about 20 degrees apart), the 
timescale distributions of microlensing candidates 
towards the two Clouds should be nearly identi- 
cal if lenses belong to the galactic halo. Also, 
the event rates should be comparable, although 
the ratio is more halo model dependent. At 
present, EROS has analysed two seasons of SMC 
data and MACHO has not yet presented its 
detection efficiency towards the SMC. From the 
published EROS efficiencies, and assuming that 
the MACHO efficiencies towards the SMC are simi- 
lar to those towards the LMC, it can be expected 
that the completed experiments will have gath- 
ered between five and ten microlenses towards the 
SMC. This should allow a significant comparison 
of the timescales (see also the discussion in |p3|). 

Finally, let us mention that, given the scarcity 
of our candidates and the possibility that some 
observed microlenses actually lie in the Magel- 
lanic Clouds, EROS is not willing at present to 
quote a non zero lower limit on the fraction of 
the Galactic halo comprised of dark compact ob- 
jects with masses up to a few solar masses. 
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